Vector Curvaton without Instabilities 
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A vector curvaton model with a Maxwell kinetic term and varying kinetic function and mass during 
inflation is studied. It is shown that, if light until the end of inflation, the vector field can generate 
statistical anisotropy in the curvature perturbation spectrum and bispectrum, with the latter being 
predominantly anisotropic. If by the end of inflation the vector field becomes heavy, then particle 
production is isotropic and the vector curvaton can alone generate the curvature perturbation. The 
model does not suffer from instabilities such as ghosts and is the only concrete model, to date, which 
can produce the curvature perturbation without direct involvement of fundamental scalar fields. 



There has been a recent outburst of activity in the 
cosmological implications of vector fields. In particular, 
the possible contribution of vector fields to the curvature 
perturbation C is receiving growing attention P, [3j 3' 
with emphasis on characteristic signatures such as sta- 
tistical anisotropy in the spectrum and bispectrum of 
C S Hi H 0]- The forthcoming observations of Planck 
satellite may well detect statistical anisotropy in the near 
future. On the theoretical side, it is a challenge to inves- 
tigate whether the curvature perturbation can be gener- 
ated without fundamental scalar fields. Even though the- 
ories beyond the standard model are abundant in scalar 
fields, no scalar field has ever been observed. The LHC 
may or may not confirm the existence of scalar fields. In 
the latter case it is imperative to find alternatives for the 
generation of (. 

The pioneering work in Ref. [l| offered such a pos- 
sibility by showing that a massive Abelian vector field 
alone can, in principle, generate C following the curvaton 
mechanism, originally introduced for scalar fields Q . The 
only requirement is that some suitable mechanism breaks 
the conformality of the vector field during inflation, such 
that an (almost) scale invariant spectrum of vector field 
perturbations is generated. This was firstly attempted 
by introducing a non-minimal coupling to gravity, of the 
form ^RA^ In another attempt, a non-trivial varia- 
tion of the kinetic function was considered 0]. However, 
in these early works on the vector curvaton, statistical 
anisotropy was ignored. But, as shown in Ref. if 
particle production of the vector field perturbations is 
strongly anisotropic the statistical anisotropy in the spec- 
trum can be excessive. In this case observational bounds 
force the vector field contribution to C to be subdomi- 
nant. This was found to be the case for the ^RA'^ model 
[1]. Moreover, the model was criticised for giving rise to 
instabilities, such as ghosts @. In view of these devel- 
opments, here we reexamine the model of Ref. Q taking 
statistical anisotropy fully into account. 
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A massive, Abelian vector field has Lagrangian density 



(1) 



where / is the kinetic function and the field strength ten- 
sor is Ff^i, = dfj,A^ — d,yA^. Eq. ([T]) can be the Lagrangian 
density of a massive gauge field. However, we need not 
restrict ourselves to gauge fields only. If no gauge symme- 
try is considered the Maxwell type kinetic term is moti- 
vated from being one of the few (three) choices [1(1 which 
avoids introducing instabilities, such as ghosts [9|]. 

We assume that the contribution of the vector field 
to the energy budget of the Universe during inflation is 
negligible, i.e. inflation is isotropic. We also assume that, 
during inflation, / = /(t) and m = m(t). 

The equations of motion for the vector field perturba- 
tions have been obtained in Ref. In Ref. 0] these 
are expressed in terms of the physical (in contrast to co- 
moving) vector field W^, whose spatial components are 
W = ^/f A/a, where a{t) is the scale factor of the FRW 
metric. To study particle production we promote the 
perturbations of the vector field to a quantum operator: 
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where k = |fc|, w\ = w\{t, k) are the mode functions and 
A = i, i?, II with L, R denoting the Left and Right trans- 
verse polarisations respectively. The polarisation vectors 
ex{k) can be chosen as 

ei = -1^(1,^,0), eR^^{l,~i,Q), e||= (0,0,1), (3) 

with k = k/k. We consider canonical quantisation with 

a\{k),a\,{k') = {2Tr)^S{k — k')Sx\' ■ The equations of 

motion for the mode functions are the same as the ones 
for the perturbations of the vector field, because the lat- 
ter are linear. By solving these equations it is found that 
all the components of the vector field perturbations ob- 
tain a scale invariant spectrum if Q 
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during inflation and ^ H, where M = m/^/f is the 
mass of the physical vector field and the star denotes the 
time when the cosmological scales exit the horizon. 
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In the following, unless otherwise stated, we focus 
on the case, when / oc a~^, which has the richest phe- 
nomenology. We will also assume that / 1 so that the 
vector field becomes canonically normalised at the end 
of inflation. Thus, if our vector field is indeed a massive 
gauge field, then the field remains weakly coupled during 
infiation because / ~ where e is the gauge coupling. 

The mode functions of the transverse components sat- 
isfy the equation Q : 



a,2 + 3ifa, + (|)'(i + r2) 



— , 



(5) 



where wl = wr = because the theory is parity invari- 
ant, and r = aM/k, with M (x during infiation. The 
solution to Eq. ^ is well approximated by 



w+ = < 



where C = 




■ifj J3ix) -iJ„^{x) 



constant and 



a; > 1 > z , 
, z > 1 > x. 



X = k/aH and z = M /3H , 



(7) 



which also implies z oc a'^ and r = 3z/a; cx a^. Subhorizon 
scales correspond to x > 1, while the field is light when 
z ^ 1 so scale invariance requires z* <C 1 . For Eq. ^ we 
used the vacuum initial condition io+ = — i=e*^. 

The mode functions of the longitudinal component sat- 
isfy the equation Q : 



Hdt 
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Using the initial condition 'W'll = -^^e", the solution 



to the above is well approximated by 
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f Ji(z)-^J i(z) 



,a;>l> z, 
x,z<^^ 
z> l>a;. 



where 7 = y 1 + 73- is the Lorentz boost factor, which 

takes us from the frame with k = (where there is no 
distinction between longitudinal and transverse compo- 
nents) to that of momentum fe 7^ 0. 

From the above we can obtain the power spectra of the 
super horizon perturbations as Vx = ^\w\\. Eqs. ^ 

and (O suggest that, when x ^ 1, there are two regimes 
depending whether the vector field is still light at the 
end of inflation (z <C 1) or not (z > 1). If the field is 
light then the typical values of the perturbations are: 
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Thus, in this case V\\ ^ V+ and particle production is 
strongly anisotropic. On the other hand, if the field be- 
comes heavy by the end of inflation, then Eqs. ([6| and 
(ini) can be written as 



2Vt iirf' ^+ = v^Jii^) + i^'Vi^-i(^) , (11) 
2\/f iirf' ^11 = (^) - Ti-^-i (^) ' (12) 

which, in the limit z 3> 1 and a; — > 0, give 

Since z ^ 1, the frequency of oscillations is very large 
compared to the expansion rate H and we can use the 
average values of the power spectra over many oscillations 
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Therefore, if the vector field becomes heavy by the end 
of infiation, then particle production is isotropic. 

As shown in Rcf. Q , the contribution to the curvature 
perturbation spectrum from a vector curvaton field is 



4^ 



V+ + (-Pji -V+)iW ■ kf 
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where W = W/W, W = \W\ and Qa 



Pa/p, where pA is the density of the vector field. As in 
the case of the scalar curvaton paradigm, Eq. (fT5|) should 
be evaluated at the time of decay of the curvaton field. 
Thus, we see that V^^ is isotropic only when particle pro- 
duction is isotropic and ~ 7^+. In the opposite case, 
can be strongly isotropic. If this is so, the vector 
curvaton can only generate a subdoniinant contribution 
to Ci with the dominant contribution coming from an 
isotropic source, such as a scalar field. The significance 
of the vector field in this case is confined to generating 
statistical anisotropy in the spectrum and bispectrum. In 
our model, the isotropy of V(^a is determined by whether 
or not the vector field is heavy by the end of infiation. 

Statistical anisotropy in the spectrum can be quanti- 
fied as [H 



Pc(k) 
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where is the isotropic part, d is the unit vector along 
a preferred direction and wc kept only the leading term 
in the anisotropy. Currently, the data provide only a 
weak upper bound on the anisotropy parameter, namely 
9 < 0.3 tni. Using the 5N formalism, we find [3 



5 = 



■Pii - V4 



(17) 



where V,j, is denotes the spectrum due to an isotropic 
source, such scalar field (j), which must be domi- 

(10) nant if 9^ 7^+. We also defined /3 = , where 
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= ^ and Na= \Na\ with = and Na 

quantify how much the elapsed e-folds are affected by the 
perturbations of the scalar and vector fields respectively. 

Statistical anisotropy in the bispectrum of C generates 
anisotropic non-Gaussianity, which can be quantified by 
the non-linearity parameter, as 



./NL = /i^^ni+5NLVl^: 



(18) 



in analogy with Eq. where = with W ± 

being the projection of W onto the plane defined by the 
three fci.2,3 vectors, which determine the bispectrum. 
For a parity invariant theory, we have 



1 + (P + 



(19) 



where k = (k = 1) in the local (equilateral) configura- 
tion and p = {V\\ - V+)/'P+. 

Let us discuss first the anisotropic case. This is re- 
alised if the vector field remains light until the end 
of inflation, i.e. m<C-ff*, where to is the final (vac- 
uum) value of the the vector field mass and de- 
notes the infiationary Hubble scale. If the field is light, 
V\\ ^ V+. Thus, the contribution of the vector field to 
C must be subdominant. This means that /3 <C 1 since 
the number of inflationary e-folds is primarily modu- 
lated by the scalar field. Hence, the anisotropy in the 
spectrum is g ~ pV^jV^ = PI z'^- For the bispectrum we 
have /jlf£ = 20^ /D, a, where we considered Q,a ^ 1 and 
took Vcji = ^ = 'P+- Since p = z^^ ^ 1 the anisotropic 
part of the bispectrum peaks in the equilateral configu- 
ration, where we find (7nl — \p^ S> 1. Therefore, non- 
Gaussianity is predominantly anisotropic in this case. 

Let us now discuss the (almost) isotropic case. This 
corresponds to rh > H^,, i.e. the vector field becomes 
heavy and begins oscillating by the end of inflation. In 
this case, V\\ ~ V+ and excessive anisotropy in V(^^ is 
avoided. Thus, we can dispense with the scalar field and 
consider = Vqj^ , i.e. the vector curvaton alone gen- 
erates C- Consequently, 0, i.e. /? » 1. Then, from 
Eqs. (fT7|) and ([Tg]) . g ~p ~ g-^i,, i.e. statistical anisotropy 
is equal in the spectrum and bispectrum and can be non- 
negligible only if TO ~ H.^, . If to ^ H^, then g = pnl = 
and C is isotropic. In this case, /nl = ./nl = 5/417^ as 
in the scalar curvaton scenario [8| . 

From Eq. (fTB|) it is evident that the magnitude of the 
contribution to curvature perturbation from the vector 
curvaton depends on the value of the zero mode W . It 
is easy to show that the equation of motion of W , under 
the conditions in Eq. (g]), is Q 



w + mw 







(20) 



during infiation, where the dot denotes time derivative. 
When / oc (i.e. M cx a^) and assuming initial energy 
equipartition between the kinetic and the mass terms in 



Eq. H]), we have 0] 

W = Wq(^^ \/2cos(z±^), (21) 

where '0' denotes initial values. Thus, the typical value 
of the zero mode scales as oc . From the above 
it can also be shown that the density of the vector field 
remains constant during inflation pA ^ MqW§ 0|. 

After inflation, / and to assume their final (vacuum) 
values and the equation of motion for W becomes . 

W + 3HW+iH + 2H^ + m^)W ^0. (22) 

From this, we can deduce that, if the field is still light 
after inflation, then pA cx: a~^, i.e. pA is diluted as radi- 
ation, in contrast to the case of a light scalar field whose 
density remains constant. When rh > Hit), the field be- 
comes heavy and oscillates with frequency to. Then, 
one can show that, on average, pA oc and P,; = 0, 
i.e. the field acts as pressureless isotropic matter and 
can dominate the Universe without generating excessive 
large-scale anisotropy, in accordance to the findings of 
Ref. [H. Thus, it can safely play the role of the curvaton. 

Working as in Ref. and assuming prompt reheating, 
we obtain [3 
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max {Fy 



min {to; H^} 



1/4 



(23) 



where fidoc is ^^^i at the time of the vector field decay, 
(^A is the curvature perturbation attributed to the vec- 
tor field, Hdom is the Hubble scale when the oscillating 
vector field dominates the Universe, if this occurs before 
decay, and Ta = h^rh is its decay rate, with h being the 
coupling to the decay products. The latter ranges as 
to/top</i<1, where the lower bound corresponds to 
decay through gravitational couplings. 

Let us consider first the (almost) isotropic case, where 
TO > H^. Then, C ~ r^docC^ as in the scalar curvaton 
scenario. This implies n^ccH^ > C,^mp. Thus, we obtain 



> 10^ GcV 



(24) 



where we used C = 5 x 10~^. Now, because M oc dur- 
ing infiation, we have to = e'^^°='=i/*, where iVosc is the 
remaining c-folds of infiation at the onset of the vector 
field oscillations. Using this and considering gravitational 
decay, it can be shown that N'^^"" = -| hiC = 4.4 0. 
Thus, the parameter space for this scenario is 



1 < m/H^ < 10^ 



(25) 



The above range is somewhat reduced if the decay of 
the curvaton is more efficient than gravitational. Indeed, 
when h ^ 1 the upper bound is reduced by 0(10^) since 
Infidec ^ 3.1, where fi^nr r*, 10~^ because a 



f miidec 0.1, wnere iidoc 
smaller ildcc violates the current observational bounds on 



osc 



NL 



5/4f7d 



Therefore, the parameter space where 
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the vector field undergoes isotropic particle production 
and can alone account for ^ is exponentially large. 

Let us consider now the anisotropic case, where 
m <^ H^. Then, it can be shown that ( ^ g~^^^^dccCA 
A lower bound on m is obtained due to the require- 
ment that the vector curvaton decays before big bang 
nucleosynthesis Ta > T^^^/mp, where Tbbn ^ 1 MeV. 
Hence, we get 

H^> ^ 10^ GeV , (26) 
while the parameter space for this scenario is 

10 TeV < m < . (27) 

Before concluding let us briefly discuss also the case 
when f oi a? in which too a scale invariant spectrum of 
all the components of the vector field perturbations is at- 
tained. In this case M = constant throughout inflation. 
Thus, we have to = M* <C if*. It turns out that Eq. (fTO|) 
is also valid in this case, but crucially 1 ^ z = constant, 
which means that 1 <C 'P\\/'P+ ~ constant. Thus, parti- 
cle production is always strongly anisotropic. The predic- 
tions for g and ^nl are identical with the case when when 
/ (X a~^. The zero mode satisfies Eq. (PH)) but initial 
energy equipartition suggests that W = Wq = constant 
during inflation, which also implies that pA ~ MqWq. 
Since the evolution after inflation is insensitive to the 
scaling of / during inflation, the bounds from the curva- 
ton mechanism are the same with those in the anisotropic 
case discussed above for / oc a"'* (e.g. Eqs. ((26)) and ([27]) 
remain valid). All in all the predictions of our model for 
/ oc a~^^^ are identical, when rh ^ H^,, while the possi- 
bility that TO > H^, can only be realised when / oc a""*. 

In conclusion, we have presented and analysed a new 
model of vector curvaton field, which docs not suffer from 
instabilities such as ghosts. The model corresponds to 
a massive Abelian vector field with a Maxwell type ki- 
netic term, whose kinetic function and mass scale during 
inflation according to Eq. The implications of our 
model depend on whether the field is light or not by 
the end of infiation If it is light, then particle produc- 
tion is anisotropic. In this case, the vector curvaton's 
contribution to the curvature perturbation must be sub- 
dominant with the dominant part due to some isotropic 
source such as a scalar field. However, the vector field can 
generate statistical anisotropy in the spectrum, which 
could be observable. It also generates anisotropic non- 
Gaussianity along a direction correlated with the statis- 
tical anisotropy in the spectrum. In fact, the anisotropy 
in the bispectrum is dominant to the isotropic part, which 
means that, if a non-Gaussian signal without angular 
modulation is indeed observed e.g. by the Planck mis- 
sion, our model will be falsified in this case. However, in 
the case when the field becomes heavy by the end of in- 
flation, particle production becomes isotropic. Then, the 
vector curvaton's contribution to the curvature pertur- 
bation can be dominant. Indeed, we can dispense with 
contributions from other sources such as scalar flelds and 



consider that the vector curvaton alone generates the cur- 
vature perturbation in the Universe. If there is some 
residual statistical anisotropy, we have shown that its 
magnitude is the same in both the spectrum and bispec- 
trum, which is a characteristic signature of this possi- 
bility. However, if the vector fleld is much heavier than 
the Hubble scale by the end of inflation, any anisotropy 
is negligible and the implications of the vector curvaton 
are indistinguishable to the scalar curvaton scenario. We 
have found that there is ample parameter space for both 
the above cases (see Eqs. (^5]) and (P7)) ). 

The curvature perturbation spectrum obtained is ex- 
actly flat. We can get a red spectrum (currently favoured 
by the data) by assuming that inflation is not exactly de 
Sitter, with e = -H/H'^ ~ IQ-^. The model needs to be 
tuned according to Eq. but this has to be contrasted 
with the unavoidable fine-tunning of scalar field models 
(77-problem) . Our model may be realised in the context of 
supergravity (supermassive gauge fields) or superstrings 
(vector fluxes in the bulk). 
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